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ABSTRACT 

In this paper we investigate the effect of a pinned superfluid component on the grav- 
itational wave emission of a steadily rotating deformed neutron star. We show that 
the superfluid pinning allows the possibility for there to be gravitational wave emis- 
sion at both the stellar spin frequency Q and its first harmonic, 2tt. This contrasts 
with the conventional case where there is no pinned superfluidity, where either only 
the 2f2 harmonic is present, or else the star undergoes precession, a feature which is 
not believed to be common in the known pulsar population. This work motivates the 
carrying out of gravitational wave searches where both the f2 and 2Q harmonics are 
searched for, even in targeted searches for waves from known pulsars which aren't ob- 
served to precess. Observation of such a two-component signal would provide evidence 
in favour of pinned superfluidity inside the star. 
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1 INTRODUCTION 

Spinning neutron stars are a source of potentially detectable gravitational waves, and so are of interest for the cur rent 
generation of laser interferometer gravitational wave detectors I Abbott et al. 2009a : Acernese et al. 20081 : Grote 20081 ) . A 
number of searches have already been Dcrformcd, so far producing upper limits on the gravitational wave emission (Abbott 
et al. l2004l . l2005al lbl. l2007a llb ll2008allb l l2009bl lci). A highlight was the result of Ubbott et alj (|2008bl) , where the energy being 
emitted in the gravitational wave channel by the Crab pulsar was shown to be no more than 6% of the total spin-down energy 
budget. 

Most gravitational wave searches have assumed emission at a single gravitational wave frequency. In the case of searches 
targeting pulsars of k nown spin frequency Q, the gravitational wave signal has been assumed to be at frequency 2fi (see e.g. 
Abbott et al. |2007ah h This is exactly the sort of emission to be expected from a steadily sp inning triaxial body, rotating 
about one of the principal axes of its moment of inertia tensor (jShapiro fe Teukolsky 1983). The necessary deformation 
could be supported by strains in the solid crust, or by strong internal magnetic fields l Ushomirskv. Cutler fc Bildsten 2000l ; 
Haskell et al.1 2008b . More complicated gravitational wave emission occurs , with multiple harmonics, if the spin and principal 
axes are misaligned ( Zimmermann fe Szedenits. Jr. 19791 : Zimmermann 1980l : Jones fe Andersson 2002 ; Van Den Broeck 
2005). However, such free precession wou l d prod uce modulations in the time of arrival, polarisation, and beam shape of the 
radio pulsations (see iJones fc Andersson! I|200lr) and references therein) . It is the absence of clear evidence of precession in 
gravitational wave candidates that has led to most searches assumin g a single monochromatic signal. A notable exception was 
the analysis for the Crab pulsar presented bv lAbbott et al.l l|2008bl ). where a small frequency band around 2Q was searched, 
motivated in part by the possibility of precession, and in part by the possibility of gravitational wave emission occurring from 
a stellar component with spin frequency slightly different from the radio-pulsation-producing outer crust. 

In this paper we revisit the issue of gravitational wave emission from rotating stars, by allowing for a superfluid component 
within the star. Ther e are strong reasons to believe that superfluids exist in neutron stars, and that they have an effect on 
the stellar dynamics (|Saulsl l 1989). Superfluids rotate by forming an array of vortices, with the number of vortices per unit 
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area being proportional to the rotation rate. An interaction between these vortices and the charged component of the star is 
believed to result in 'pinning' of these vortices, preventing that portion of the superfluid that is pinned from participating in 
the smooth spin-down of the star. Instead, the spin rate of the pinned superfluid remains constant, so that a spin lag builds 
up between the superfluid and the pinned vortices. When this lag becomes sufficiently large, some mechanism is believed to 
result in a catastrophic unpinning, spinning up the star's charged component, and producing one of the well known pulsar 
glitches. 

As was pointed out bv lShahaml (|1977l ). in addition to producing glitches, a pinned superfluid component also has important 
repercussions for the precessional dynamics of a star. Classically, the free precession of a rigid body is determined by the 
principal moments of inertia I2, H)- However, Shaham showed that a pinned superfluid alters the motion, effectively 
acting as a large gyroscope, 'sewn' into the body of the star. If the star is set into precession, the angular momentum vector 
of the pinned superfluid follows the wobbling motion of the star, and this must be included in the solution of the Euler rigid 
body equations of motion. 

The detailed form of the free precession that then occurs is given in IShahaml l| 19771) . We will not concern ourselves with 
the corresponding gravitational wave emission here. Instead, we will confine our attention to the allowed non-precessional 
motion of such a star, and the consequent gravitational wave emission. The key point is this: if one allows for the possibility 
of the superfluid pinning along an axis that isn't one of the principal axes of the star's moment of inertia tensor, even the 
non-precessing motion is interesting: despite rotating steadily at rate fi, there is gravitational wave emission at both f2 and 
The relative strengths are dependent on the moments of inertia and relative orientations of the stellar components. Basically, 
this motion is possible because the moment of inertia of the pinned superfluid is likely to be many orders of magnitude 
greater than the asymmetries (I2 — I\,Iz — 12,1s, — h), which (as we show below) means that, to a good approximation, 
the non-precessing motion corresponds to rotation about the superfluid pinning axis, not to rotation about any one of the 
principal axes. This means that even in the absence of any observed precession in the radio data, a given star might possibly 
be radiating at both Q. and 2Q, motivating the carrying out of gravitational wave searches that account for thi s . 

There are two different locations where superfluid pinning may occur in a neutron star (see e.g. Saulsl 1 1989h ). One 
location is the inner crust, the other is the core, and the physical conditions in the two reg imes are very different. In the crust, 
pinning may occur because of an interaction between the vortices and the nuclear lattice IjAnderson fc Itohlll975T ). In the case 
of small relative velocities between superfluid and vortex, the pinning may be strong enou gh to prevent motion of the vortices 
with respect to the lattice; an explicit model of this was presented recently by Link! ( 2009T ) . This crustal superfluid is believed 
to make up only a few percent of the star's total moment of inertia. The other possible pinning location i s the c ore, where 
there is believed to be a strong interaction between the vortices and magnetic flux tubes ( Ruderman 19761 . 1991al lbh. In this 
case, the bulk of the star's moment of inertia might be a pinned superfluid. For the most part, the exact location and cause of 
the pinning is not important for our purposes — what matters is that the pinning exists and is maintained for periods longer 
than the gravitatio nal wave observation time. Th e long separation between glitches in young pulsars indicates that this may 
indeed be the case i Lvne fc Graham-Smith 19981 ). 

The plan of this paper is as follows. In section [2] we present our basic rigid body stellar model and derive the form of 
the non-precessional solution. In section [3] we calculate the corresponding gravitational wave emission, count the number of 
parameters required to describe the signal, and discuss the angular pattern of the emission. In section [4] we show that the 
results obtained still apply when the assumption of rigidity is relaxed to one of elasticity. In section [5] we discuss what might 
limit the wave emission and what mechanisms might prevent our model being realised. In section [6] we talk about two other 
(non-superfluid) mechanisms that can lead to radiation at multiple frequencies, and show they are likely to be less important 
than the superfluid pinning mechanism. Finally, in section [7] we speculate as to how the sorts of deformed stars considered 
in this paper might arise in nature. 



2 BASIC MODEL 



We will begin by following IShahaml (|1977f ). extending his treatment slightly by modelling a star consisting of a triaxial rigid 
'crust' and a pinned superfluid component of spherical moment of inertia (Shaham assumed a biaxial crust). Note that here, 
and for the entirety of this paper, by 'crust' we mean not just the outer solid phase but all parts of the star apart from 
the pinned superfluid. Referring our equations to the crust's body frame and denoting by {n^, n^, n^} unit vectors along the 
crust's principal axes, the crust's moment of inertia tensor is 



r C rC 1 1 . r c 2 2 . r c 3 3 
hb = h n a n b + I 2 n a n b + I 3 n a n b , 

while the superfluid's is 

t-SF r SF s 
lab ~ 1 a b- 

The star's total angular momentum is then the sum of the angular momenta of the two components: 



(1) 



(2) 
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j a =ignf+i s ^, (3) 

where Q.^ and fif F are the crustal and superfluid angular velocities, respectively. Writing out the components of J a explicitly: 

Ji = + i SF nf F , (4) 

J 2 = 7 2 c fi 2 3 + 7 SF fif, (5) 
J 3 = 7 3 c ^ + 7 SF r!f. (6) 

If the pinned superfluid has a rotation axis given by the spherical polar angles (9 SF , SF ) with respect to the crust frame then 

Of = fi SF (sin0 SF cos0 SF ni + sin SF sin S V ,+ cos 9 SF n 3 a ). (7) 

The Euler equation of motion is 

<IJa +e abc nUc = 0, (8) 



dt 

where the time derivative is evaluated in the crust frame. Writing out the components explicitly we obtain: 

+ n%n%(is - 1?) + {ri2 cos6 SF -n.f sme SF sm(j) SF )i SF n. SF = o, (9) 

I2Q2 +^3^i(I? -l3) + (^3sme SF cos(j} SF -Q^cose SF )I SF Q. SF = 0, (10) 

+nfn 2 : (7 2 c -7 1 c ) + (^sin6l SF sin^ SF -r2 2 c sin6l SF cos0 SF )7 SF fi SF = 0. (11) 

The non-precessing solution is obtained by setting dj a /dt = in equation ©, which in turn implies e a (, c f2^ J c = 0, i.e. is 
parallel to J a , so that the crust rotates steadily about the fixed angular momentum vector. Equations (l9|)- (llip then reduce to 

n2n3(7 3 3 -7 2 c ) + (fi 2 3 cos6l SF -n3 C sin6l SF sin^ SF )7 SF n SF = 0, (12) 

n3 J nf(7 1 c -7 3 c ) + (^sine SF cos^ SF ^^cose SF )7 SF n SF = 0, (13) 

fi?fi 2 (7f - 7?) + («?sin0 SF sin<£ SF - f2 2 c sin SF cos SF )7 SF fi SF = 0. (14) 

If we reg ard fi SF ,6» SF ,<?i> SF as fixed, as is appropriate for a perfectly pinned superfluid, these equations then have a one- 
parameter family of solutions, which we will parameterise in terms of Q.3 . Then we find 

oC oC 7 SF » SF sinfl SF cos^ F 

_ S S3 in C +/ SF n SF cos 0SF> ^ 

a _ p c 7^ F ^ F sin^ F sin^ F 
having defined 

Jsi = Is-I?, (17) 

7 32 = 7 3 c -7 2 c , (18) 

7 21 = 7 2 c -7 1 c . (19) 

To fix the orientation of the body with respect to the inertial frame, define a set of three Euler angl e s (0, < t>, ip) which connect 
the body frame to the inertial frame in the standard manner, as described in Landau fc Lifshita (197a) and illustrated in 
Figure □ Then 

Qi = ^ sin 6» sin + cost/>, (20) 

2 = </> cos # cos^/> — dsmi/j, (21) 

^3 = <j>cos8 + ip. (22) 

Our non-precessing solution has all components of Qa constant, and so equations (|20|l and (|21[) imply that 9 and ip must be 
constant, leaving 

fii = 4> sin6» sin ip, (23) 

fia = 0cos<9cosi/>, (24) 

0,3 = <p cos 9. (25) 

These can be inverted easily to give 

tantf = [( l c J , (26) 
i.l 3 

4> = n c t + cp , (27) 
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Figure 1. The orientation of our body is specified by the three standard Euler angles (0, <f>, ijf), as labelled above. The fixed inertial-frame 
axes are denoted by (x, y, z), while the rotating body-frame axes are denoted by (X, Y, Z). The so-called line of nodes, N, lies along the 
intersection of the xy and XY planes. 



tantp 



na- 



ps) 



The angles 9 and i[> are constant, while is a linearly increasing function of time — it is this angle that generates the rotation. 
Elimination of nf and 0% in favour of n^ using equations (|15[1 and (|16[1 gives 



tan t 



tan t 



flit < tan 



/ 3 in c 



+ sin 



1 + 



/3 2 ng 
/ SF nf F 



1/2 



i + 



/ SF nf 



• 2 ,SF 

+ sin (f> 



h2^ 

/ SF nf F 



1/2 



+ 1 



+ <f>0, 



taxxip = tan(7r/2 — <^> SF ) 



1 + 



i SF nf 



i + 



/ SF nf 



(29) 



(30) 



(31) 



where 4>o is a constant. Clearly, there is a rather complicated relationship between the orientation of the body with respect 
to the inertial frame and the parameters Q§ , Q SF , 9 SF , (f> SF ■ 

However, we can simplify the above. The asymmetri es /31 , J32 are s ource d by crustal strains or magnetic fields, and so 
are likelvtc^jDe^er^ i sinalL > Tr^tlie i ^ strains, Haskell et al. (2006), following the^orl^o^Ushomirsky, Cutler & 

Bildsten lj2*00*oi r'*e'stonlrte'd^ 



AI 



3 x 10" 



(^brcak \ 
0.1 J 



TO 45 gem 2 ■-■"*— V 0.1 J ' ^ 

where iibroak denotes the breaking strain of the crust. Recent molecular dynamics calculations ( Horowitz fc Kadau 2009h 
estimate this breaking to strain to be approximately 0.1. Note that the above crustal deformation is an upper boun d as it 
assum es the crust is maximally strained. In the case of an internal magnetic field of strength ~ B a naive estimate is ( Jonesl 
20021) 



AI 



10 45 g cm 2 



10" 



uo 12 cy ' 



(33) 



although more detailed calculations show the result depend s on the assumed field geometry (jBonazzola fc Gourgoulhonl 1996; 
Haskell et aflbood Iciolfi et alj|2009l ; lLander fc Jonesl 2009) and on the superconducting nature of the core (see lCutlerl (|2002l ) 



and references therein). 

On the other hand, the m oment of inertia of the superfluid is relatively large: Isf/I 
for core pinning (|Saulslll98Sh . We can therefore make use of the inequalities 



10 2 for crustal pinning, I SF /I ~ 1 



© 2009 RAS, MNRAS Q00.[THT9l 



GW emission from superfluid neutron stars 5 



ISFQSF « 7 SF n SF « 

to simplify our expressions for the Euler angles 

6» w 8 SF , (35) 
cos 8 bt 

iP « ^/2-0 SF , (37) 

This shows that the star rotates about an axis very close to n„ F . So, after having gone through a careful analysis, the basic 
picture is very simple: the star rotates steadily about an axis very close to the superfluid pinning axis. From this point on we 
will simplify our notation slightly, writing (fi„ fi„ ) 1//2 simply as f2, the observed stellar spin frequency. 



3 GRAVITATIONAL WAVE EMISSION 

3.1 Calculation of the gravitational wave signal 

Having obtained the Euler angles (8, cp, ip) giving the orientation of the body with respect to the inertial frame it is straightfor- 
ward to calcula te the correspo nding gravitational wave emission. To do so we will calculate the source's multipole moments, 
as described by Thornei 1 198ut) ; this will make clear the separation between the frequency harmonics and will also be useful 



in discussing maximum possible wave amplitudes in section 13.31 The multipole moments of the source are complex scalars 
denned by integrals over the star's density field p: 

Qim = J pr l Y* m d 3 x. (38) 
When the density field is written as the sum p = Y>8pi m Yim this reduces to 

Qlm = J Spi m (r)r l+2 dr. (39) 
When the body has orientation 8, cp, ip relative to the inertial axes we find 

1 / 15 _i0r 



Qui = -^—e-" p [-l2ism2^sm9-ism29(I 2 icos^^h 1 )], (40) 
Q22 = — -1 /^— e _2l ^[/2i(cos 2 ip cos 2 8 — sin 2 ip) + J31 sin 2 9 — H21 sin 2ip cos 01. (41) 

2 y 27T 

The Q21 multipole moment generates radiation at the spin frequency Q, while the Q22 multipole moment generates radiation 
at 2Q. 

In calculating the gravitational wave emission, we will consider an observer a distance r from the star whose coordinate 
system (X, Y, Z) is obtained from the inertial (x, y, z) frame by a rotation through an inclination angle l about the inertial 
a;-axis, followed by a translation through a distance r along OZ. The gravitational wave then propagates along OZ and, within 
the mass quadrupole approximation to general relativity, the two polarisations which make up the transverse traceless wave 
field are 

h + = sin l cos t{ [I21 cos 2 ip — sin 28 cos <j) — I21 sin 2tp sin 8 sin 0} , (42) 

h x — sin t{ I21 sin 2tp sin 8 cos 4> + (J21 cos 2 ip — /31) sin 28 sin <fi} , (43) 

h + — (1 + cos 2 t){ [J21 (cos 2 ip cos 2 8 — sin 2 ip) + J31 sin 2 9] cos 2<f> — I21 sin 2tp cos 8 sin 20} , (44) 

h x — 2 cos ij/21 sin 2tp cos 8 cos 2<p + [(I21 (cos 2 ip cos 2 8 — sin 2 ip) + I31 sin 2 8] sin 2<ji} . (45) 

These equations, together with the equations (|29[l - (|31[l ( or alternatively the approximate equations (|35|l - (|37p l then give a 
complete specification of the gravitational field. As anticipated, there are harmonics at both SI and 2Q. 

3.2 Counting the number of parameters 

The actual signal received by a detector can be written as 

h(t) = F+(t)h+(t) + F x (t)h x (t), (46) 
© 2009 RAS, MNRAS 000.HM91 
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where F+ and F x are the beam pattern factors whose form can be found in Jaranowski, Krolak fc Schutz (1998). They depend 



upon the detector location, source position on the sky, and on the orientation of the source's spin axis. This introduces three 
unknown parameters over and above those appearing in the expressions for h+ and h x , specifically the source's location on 
the sky (which can be specified by right ascension a and declination 5), and an angle giving the projection of the source's 
spin axis on the plane perpendicular to the line-of-sight (which can be specified by the polarisation angle Vw)- 

We can now count the number of parameters required to reconstruct the full signal. Let us first examine two special cases, 
to demonstrate consistency with previously published results. 

If we put 9 = in (|42)l - (f45)) we obtain 

2Q, 2 

h+ = — I 2 i(l + cos 2 i) cos 2(<t> + ij>), (47) 

2Q 2 

h x = J 2 i 2 cos i sin 2(<^) + ?/>). (48) 

This corresponds to a triaxial star spinning abou t a prin cipal axis; it is the sort of signal that has been assumed in many 



gravitational wave searches (see e.g. lAbbott et al.l (12007af0 . A possible choice of parameters to describe the signal is: 



2 hi 



-,Q,((f> + il))o,L, ippoi, a,6 > . (49) 



In this case there are 7 parameters (note that in this case the angles <f) and ip are degenerate, and only their sum appears in 
the waveform). There is radiation only at the 20. harmonic. 
If we put 7 2 i = in (|42)) - (f4"5)) we obtain 



ft+ = — sin l cos 1/31 sin 29 cos cj>, (50) 



Q 



= — — sin t/3i sin 20 sin 0, (51) 
r 

h 2 ^ 1 — (1 + cos 2 t)/3i sin 2 9 cos 2(j>, (52) 

r 

h 2 ^ 1 = 2 cos 1/31 sin 2 9 sin 2(f>. (53) 

r 

This is identical to the wave field of a processing biaxial star, and therefore is of exactly the same form as given in Zimmermann 
& Szedenits, Jr. 1 19791 ) and Jaranowski. Krolak fe SchutzT i 19981) . Note, however, that in our case the motion is not one of 







free precession — the star rotates steadily about the fixed angular momentum vector, with no superimposed rotation about 
the crust's symmetry axis. A possible choice of parameters is: 

<f> o ,0,i,ip po i,a,6 > . (54) 

r J 

In this case there are 8 parameters, with the angle 9 being the extra parameter as compared to the case considered above. 
There is radiation at both the O and 20. harmonics, and 9 controls the relative strength of these. In the limit of small 9, they 
scale as 9 and 9 2 respectively, so the O harmonic can dominate, but both go to zero (at fixed /31) in the limit — > 0. 

Finally, consider the general case. The wave field is then as given in equations (|42|| - (|45|l . A possible choice of parameters 

is 

O 2 121 I31 ] 

, 1—!^, 0o, #,■*/'>''> i/j po i,a,5 > . (55) 

r hi J 

In this case there are 10 parameters, with the ratio /31//21 and the angle if) being the extra parameters as compared to the 
case above. In general, there is radiation at both fi and 20, but the relative strengths depends on the parameters in a rather 
complicated way. Let us consider the extremes. There are solutions where the star radiates only at 2Q, with the signal being 
proportional to one of iai, /31, J32. These correspond to the familiar case of rotation of a triaxial star about a principal axis. 
More interestingly, there also exist solutions where the star radiates only at Q. The solutions are: 



>.n> <■ = (). taiifl = IQ21I = vS V_/l ' :7 - :i - 



hi , /15 



cosi/> = 0, sm9 = x /+^-, |Q 21 | = J (57) 



cos# = 0, sinV> = ^f^' IQ21I = y ^vTaTfe. (58) 

Clearly, to obtain radiation only at the Q, harmonic, special values of the parameters are necessary. For general values of the 
parameters, there will be radiation at both harmonics. 
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3.3 Gravitational wave amplitudes 



Let us now examine the relative strengths at which the two quadrupolc moments can radiate. As noted bv Ushomirsky, Cutler 
& Bildsten l^OOoh . the Q22 perturbations are somewhat more efficient gravitational wave emitters, in the sense that for a 
given magnitude of the quadrupole moment (i.e. IQ21I = IQ22I), the energy flux in the 2Q radiation is greater than that in 
the Q radiation. To see this we can calculate the inclination- angle (t) dependent gravitational wave energy flux using (|Thorne! 



198C 



F = 



32tt 



;tt;tt 



(59) 



where the angle brackets denote a time average over several periods. This gives the following when written in terms of multipole 
moments: 



= 

F 2n (0 = 



1 Q b 



-^-T\Q2i\ 2 sin 2 i(l + cos 2 t), 
4 K I O22 1 2 [( 1 + cos 2 1) 2 + 4 cos 2 l] . 



15 r 2 

When integrated over the sphere the corresponding total luminosities are 

47T 

75" 
256tt, 



En 
E2C1 



-o 6 |Q2i| 2 , 



75 



-n 6 |g 22 | 2 . 



(60) 
(61) 

(62) 
(63) 



As noted bv lUshomirskv. Cutler fc Bildsten (2000), for a given magnitude of the quadrupole moment, a Q22 perturbation is 
64 times more luminous than a Q21 perturbation. 

However, in terms of detectabilit y, the relevant quantity is not the ene rgy flux but the signal-to-noise ratio, which is linear 
in the gravitational wave a mplitude ( Jaranowski, Krolak fc Schutz 1998[ ). For a star radiating continuously at gravitational 
wave frequencies Q, and 2f2, Jaranow'skL*l^r6l£dTfc*ScnutJ i 19981 ) showed that the signal-to- noise is 



d = y/dj + d 
where 

9 

dx = 



2 



hl(t)dt, 



h 2 (t)dt. 



(64) 

(65) 
(66) 



In this equation, the time integral is over the duration of the observation, Sh{f) is the detector's noise spectral density, and 
h\(t),h2(i) are defined in terms of the beam functions of equation (|46[) : 



hi(t) = F+(t)h i i(t) + F x (t)h l i(t) 
h 2 (t) = F+(t)hf(t) + F x (t)hT(t). 



(67) 
(68) 



The wave amplitudes ft? (i) , h% (t) , hf (t) , hf 1 (t) are given in equations (I42|) - (I45|I . In the case of gravitational wave searches 
at a single frequency from a known pulsar, it was found that a signa l-to-noise threshold of 11.4 was necessary to achieve a 
false alarm probability of 1% and a false dismissal probability of 10% (| Abbott et al.ll2004il . This signal-to-noise threshold will 
have to be recalculated for a two-frequency search; we will discuss this further in section [7] 

The value of the signal-to-noise for a particular star, and the relative importance of the two harmonics, depends upon all 
10 of the parameters given in equation (155[) . In order to gain a rough idea of the relative importance of the two harmonics, 
we will define an approximate signal-to-noise ratio for each according to 



P = "v j °"° _ ( 69 ) 

\J <Sfe(/Gw) 

where T b s is the duration of the observation (see I Jaranowski. Krolak fc Schutz (1998) for a justification of the scaling with 
observation time) . For the gravitational wave amplitude h we will take the square root of the time-averaged sum of the squares 
of the plus and cross gravitational wave polarisations, i.e. 

h=< (h + ) 2 + (h x f > 1/2 , 

where the angle brackets denote time averaging. Then equation (|59|) gives 

4 



(70) 



h=< (h+) 2 + (h x ) 2 > 1/2 = 



(71) 
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Figure 2. Contour plot showing normalised signal strengths of the Q (inner contour) and 2Q (outer contour) harmonics, assuming 
IQ21I = I Q22 1 ■ For the H radiation we plot sint(l + cos 1 ) 1 / 2 , for the 2Q radiation 2[(1 + cos 1 ) 2 + 4cos 2 t] 1 / 2 . The rotation axis (t = 0) 
lies along the vertical. 



allowing us to write the signal-to-noise ratios for the two harmonics as 
A 

2A 



pn(i) = A J Q21 sm t (l + cosQ 1/2 , (72) 
2 A 

pan(i) = lQ 22 |[(l + cosQ 2 +4cos 2 t] 1/2 , (73) 



where 

, 1 / 2 



A = 2 



(~) Tj/^ 2 . (74) 

To compare the relative signal strengths of these harmonics, the inclination angle-dependent coefficients of A\Q\ 2 / Sh(f) 
appearing in the above equations are presented in Figure [2] The Q21 multipole emits preferentially in the rotational equator 
(t = t/2), emitting no radiation along the rotation axis (t = 0), while the Q22 emits preferentially along the rotation axis, 
with a non-zero wave amplitude in all directions. Clearly, for all observer locations and for [Q21I = IQ22I, the 20. radiation 
is stronger than the O radiation. The ratio of the sky-averaged gravitational waves strengths (as denned in equation (|71[) ) is 
approximately 3.9 ~ 4. Of course, the ratio of the possible maximum signal-to-noise ratios would also fold in the frequency 
dependence of the detector noise between frequencies O and 20.. 

If instead we were to place an upper limit on the wave amplitude of a star known to be losing kinetic energy at a rate E, 
then equations (|62 ll and (1631) show that the upper limits on the multipolcs stand in the ratio IQ™ ax | = 8]Q?9 ax i (Ushomirsky, 
Cutler & Bildsten boOo T It then follows that a star losing energy at a known rate is, when sky-averaged, 8/3.9 w 2 times 
stronger in gravitational wave amplitude if it emitting via a pure Q21 multipole than when emitting via a Q22 multipole. Such 
an upper limit is appropriate for known pulsars of measured spin-down rate (if one assumes all of the kinetic energy being 
lost is converted into gravitational wave emission) , or for low-mass X-ray binaries of known X-ray luminosity (if one assumes 
the spin-up accretion torque, as estimated from the X-ray luminosity, is balanced by a spin-down gravitational torque). 



4 EXTENSION TO AN ELASTIC STAR 

It is straightforward to allow for the rotationally-induced deformation of the star by adding terms to the moment of inertia 
tensor. For the crust we amend equation to 

/C rC_ 1 1 1 t-C 2 2 1 tC 3 3 i a tC C C /tct\ 
ab = h n a n b + 1 2 n a n b + 1 3 n a n b + Al n n a n b , (75) 

where is a unit vector pointing along S7„ . Physically, the term AlQn^n b represents a deformation biaxial about Q.^, and 
A/q is numerically equal to the difference in moments of inertia about axes parallel and orthogonal to the f2„ induced by 
the rotation. The terms in if , if an d I3 give the moment of inertia of the star in the absence of rotation. The differences 
^31, ^32, 121 are defined as before (equations (|17|) - (|19|) ) and can be non-zero only by virtue of non-rotational forces, e.g. elastic 
or magnetic strains. 

A similarly decomposition applies for the superfluid: 
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7-SF 7-SFr i a rSF SF SF /r, fi \ 

hb = I o ab + AI n n a n b , (76) 
where n| F is a unit vector along Slf F . The star's angular momentum components then become 

Ji = (i? + MR)n? + (i SF + M s n F )^ F , (77) 

J 2 = (I 2 c + A/£)^ + (/ SF + A/£ F )nf , (78) 
J 3 = (I? + Al£){l% + (I SF + AI^Qf . (79) 

These are identical in form to the corresponding rigid star equations, (|4])-([6]), the mapping from the rigid star to elastic star 
equations being — > I± + AIq (and similarly for if, if) and I SF — > I SF + A7fj F . It follows at once that the solutions, both 
precessional and non-precessional, are of the same form for this more realistic elastic model as for the rigid star. 

There will however be a change in the corresponding gravitational wave emission, as we must now include the contributions 
from the mass asymmetries associated with the AIq and A7q F components. For the non-precessional solution it is easy to 
convince oneself that the impact on the gravitational wave emission is minimal. The fl^ vector points along the invariant 
angular momentum axis, so there is no gravitational wave emission from the A/g term. The A7|j F piece represents a biaxial 
deformation about fif F , which is misaligned with the angular momentum axis. This piece therefore does radiate. Rather than 
calculating this < 
in a cone of half- 
algebra leads to 



calculating this contribution fully, we can make a simple estimate. This deformation is of size AI^ F and moves at a rate Q, 
in a cone of half-angle cos -1 (nin a F ) about n J a , the unit vector in the direction of J a - Using the results of section [2] a little 



_. , SF sin0 SF [cos 2 SF (/3icos 2 SF + /3 2 sin 2 ^ SF ) 2 + (7 2 isine SF sin0 SF cos^ SF ) 2 ] 1/ ^ 
cos (n a n a ) » JSf^sf ' ( 80 ) 



where we have worked to leading order in the small quantities of equation (|34|) . We see that Q a F lies very close to J a ; for 
Q§ w nf F we obtain 

-i/ j sfn A/ 

cos (n a n a ) ~ w, (81) 
where we have defined 

A/ = max(/ 2 i,/3i,/32). (82) 

It follows that the A7q F term produces a gravitational wave field of a similar form to the /31 term in equations (|50[) - (153p . but 
with the replacements 731 — ► A/q F and 9 — > A7/7 SF . It follows that the contribution of this component to the gravitational 
wave emission is weaker than that of the crust as calculated in equations (|50|I - H53|) by a factor AJq F // sf . For all realistic 
scenarios this ratio will be much less than unity, so we will therefore not consider this component of the wave emission in the 
remainder of this paper. We therefore see that allowance for rotational deformation of the non-precessing star has minimal 
impact on the gravitational wave emission. 



5 LIMITATIONS AND CAVEATS 



What limits the maximum gravitational wave amp l itude ? In the absence of the pinned superfluid component the answer has 



been obtained by Ushomirskv. Cutler fc Bildstenl 1 2000l) : the Q22 multipole is limited by the breaking strain iibrcak of the 
crust: 

2 ( ^brcak 



Q 



1.2 x l(r M gem' 



V 10- 



(83) 



(see also Irlaskell et all (|2006|) ). Also, as discussed in the conclusions section of I Ushomirskv. Cutler fc Bildstenl l|2000l) . an 
identical limit is placed on Q 2 i. (That these bounds are identical can be understood by noting that the real parts of Y22 
and Y21 density perturbations are identical up to a rotation of axes). So, if the ability of the crustal strain to produce the 
asymmetries 7 2 i, /31, -Z32 were the only factor, the upper bounds on the signal strengths would be given by equation (|72p with 
IQ21I — Qmax for a pure Q 2 i deformation, and by equation (|73|l with |Q 22 j = Q max for a pure Q 22 deformation. In the general 
case of a mixed Q21, Q22 deformatio n, both perturbations contribute to the strain field, and the limit on each is reduced 
1 Ushomirskv. Cutler fc Bildsterjlioooh . 

However, the presence of the pinned superfluid complicates this picture. In our model, there is an angular velocity 
misalignment between the superfluid and the rest of the star, as calculated above (equation I81[). It follows that over the 
pinning region the crust exerts a force on the superfluid, supplying the torque required to allow for its time- varying angular 
velocity. As is demonstrated in appendix [A] this torque is supplied by the Magnus force, generated by the flow of superfluid 
past the pinned vortices. It follows that there ar e new failure mechanisms as compar ed to the c ase studied bv Ushomirsky, 
Cutler & Bildsten (2000). Firstly, the pinning force interaction between the solid crust and the vortices has to be sufficiently 
strong that the vortices don't unpin. Secondly, the corresponding equal but opposite force on the crust will strain it further. 
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This additional strain might push the crust closer to, or beyond, its breaking strain. This will also produce some change in 
the moment of inertia tensor, potentially modifying the gravitational wave emission. Finally, there is also the possibility that 
relative flow between the superfiuid and the vortices may trigger an instability. Let us make some simple estimates of the 
importance of these effects. 

To do so, note that we could set our star up with an arbitrarily large difference between S7„ F an d ■ However, in the 
physically most plausible case, these two vectors will be of equal magnitudes, but misaligned by the angle cos -1 (nfnjp), 
producing a vortex-superfluid relative velocity of order 

Av ~ MR cos _1 (n^ p ), (84) 



where R is the stellar radius. Equation (|8f I) then leads to 

Av-QR^L. (85) 



A„ ~ 6 x 10 W 1 i^L\ 



Parameterising: 

where / = fl/(2ir), and I denotes the total stellar moment of inertia. This will be useful in interpret ing some of the estim ates 
that follow. This can equivalently be restated in terms of a pinning force per unit length of vortex ( Link fc Cutler 2002): 

t SF SF/ SF C\ / Q ~\ 

f a = p K.e abc n b [v c - v c ) (87) 
where p SF is the superfiuid density and k the quantum of vorticity. Writing the relative velocity using equation (|85|) gives 



g cm 3 J \ 10 6 cm 



R 



(88) 



These equations give the relative vortex-superfluid velocities, or, equivalently, the pinning forces per unit length of vortex, 
required to sustain the non-precessional solutions described in this paper. 



5.1 Breaking the pinning 

First let us ask if the pinning is strong enough to withstand the Magnus force. There exist both observational and theoretical 
inputs to this issue. On the observational side, glitches and neutron star precession may tell us something about the pinning 
strength, while, on the theoretical side, there has been some microphysical modelling of how vortices interact with the rest of 
the star. 

Let us start with the observations of glitching pulsars. The smaller glitches, such as those seen in the Crab, may be caused 
by 'starquakes', i.e. changes in the elastic strain of the crust, and so may not carry information about the superfiuid pinning. 
However, the larger so-called 'giant glitches', such as those in the Vela pulsar, are believed to represent sudden unpinning events. 
An estimate of the relative superfluid-vortex velocity at unpinning comes from assuming angular momentum conservation 
( Link fc Cutierll2002l : Iciampedakis fc Anderssor]2 009) . ■ f the glitch represents a sudden angular momentum exchange between 



the pinned superfiuid and the rest of the star (in our language the 'crust'), then the changes in angular frequency of these 
two components are related via J SF Af2 SF + I c A£l c = . Observations indicate that, when averaged over many glitches, a 



few percent of the total spin-down is reversed in glitches ( Lvne fc Graham-Smith 19981) , corresponding to I SF /I c ~ 10 2 . It 
follows we can estimate the pre-glitch velocity lag by 



rC 

Av w AQ, SF R w — AQ, C R, (89) 

ISF 

where the quantity AQ. C is obtained directly from glitch observations. Parameterising with the Vela in mind: 



/ 1Q- 2 \ / AO c 
\I SF /I J \7 x 10" 5 Hz 



A^ 7 xl0 3 cms-M^l (-——-). (90) 



Comparing with equation (|86[) . we see that this is smaller than the relative velocity needed to sustain our non-precessing stars, 
for at least some part of the potential parameter space (viz rapidly spinning stars with large deformations and low pinned 
superfiuid fractions). However, it is likely that (|90[) is too conservative an estimate on the critical velocity for unpinning. 
Firstly, the estimate assumes that the velocity difference is reset to zero at the glitch. If only some of the velocity difference is 
relieved, a larger critical velocity would be obtained. Secondly, glitches presumably involve the most weakly pinned superfiuid 
(in a spinning down star, the weaker pinning sites will break before the stronger ones), whereas it is the strongest pinning 
locations that matter in building our non-precessional solutions. So, it is not clear if the relative velocity of equation (190ft 
limits the allowed vortex-superfluid velocities in neutron stars. 
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Turn now to free precession. The free precession per iod of a neutron star in the absence of superfluid pinning is given 
approxi mately by Pf a ~ PI/AI ( Landau fc Lifshitz 1976b . while the precession period with pinning is approximately Pf p ~ 



PI /I ( Shahamlll977l ). Given that we expect A7 <C / , it follows that the precession period is much shorter in the case 



of superfluid pinning than without pinning. If follows that an observation of free precession would yield an insight into the 
stellar interior. Unfortunately, free precession does not seem to be a common pheno menon in the pulsa r population. The 
most convincing observation is of approximately 500 day precession in PSR B1828-11 (jstairs et alJl200Ch . but more recent 



observations of this and other objects may cast doubt on the precessional interpretation (M. Kramer, personal communication). 
Nevertheless, if interpreted as precession, the lo ng 500 day perio d implies that there is no significant pinned superfluid 
component in the precessing star. As discusse d by Link (2003, 2006 ), this may h a ve significan t repercussion for the nature of 
superconductivity in the stellar core (but see Glampedakis. Andersson fc Jonesl 1 20081 2009h for a counterargument). What 
does seem likely is that, if PSR B1828- 11 really is precessin g, the precessional motion itself must have set up a velocity field 
that has caused large scale unpinning (|Link fc Cutler! I2002T ). It follows that one can estimate the critical velocity to up-pin 
even the strongest pinne d vortices by examining the relative flo w in a precessing neutron star with superfluid pinning. The 
velocity field is of order ( Glampedakis. Andersson fc Jonesl 2009h 



Av < QR6- 



(91) 



where 9, the 'wobble angle', can be estimated from observations (|Jones fc Anderssonll200lt iLink fc Epsteinll200ll) . Parame- 
terising in terms of PSR B1828-11: 



Au(PSRB1828 - 11, precession) < 8.1 x 10 3 cms 1 ( 



R 



10 6 cm/ V3° 



I SF /I 
lO" 2 



(92) 



Comparing with equation (|86|l we see that this critical velocity is low enough to be problematic for some portion of the 
parameter space of interest for our non-precessing solutions, specifically very rapidly rotating stars with small superfluid 
fractions. 



Finally let us turn to theoretical estimates of the strength of pinning. In the case of pinning in the inner crust, ILink 



( 2009h recently examined the problem of how a vortex embedded in a nuclear lattice responds to a superfluid flow. Link 
solved the equation of motion of a single vortex, including tension, Magnus and drag terms. He found that when the relative 
vortex-superfluid velocity lay below a critical value the vortex essentially remained bound to the nuclei, i.e. the perfect pinning 
regime applied. Using recent calculations of the vortex- nucleus interaction energy ( Donati fc Pizzocherol 20061) Link estimated 

Aw(crust, theory) « 10 6 -10 7 cms" 1 . (93) 

Comparing (193[) and (|86[) indicates that, by Link' s calculation at least, our non-precessing stars are safely in the regime of 
perfect pinning in the inner crust. As discussed by Link (2009), this velocity difference is large enough to allow the necessary 
accumulation of differential rotation to explai n giant glitc hes, but too large to be consistent with an observation of slow free 
precession. In the case of pinning in the core. iLinkl |2003J) estimates a maximum pinning force of 3 x 10 15 B\!^ dyn cm" 1 , or 
in terms of a relative velocity, 

At;(core, theory) « 5 x 10 3 cms" 1 f — 5— ) ^ ( 3 x 1Ql4 g cm ) . (94) 



Comparing with equation (|86jl we see that, by this estimate at least, the core vortex-fiuxtube interaction may not be strong 
enough to sustain pinning for non-precessing stars with high spin frequencies and small superfluid fractions, and the problem 
is more severe for weakly magnetised stars. 

To sum up, the critical unpinning velocities implied by observations suggest that unpinning may be a problem for our 
non- precessional solutions, but only for very rapidly spinning stars with relatively small pinned superfluid fractions. Slower 
spinning stars, or stars with larger pinned superfluid fractions, should be immune from Magnus-force generated unpinning. 
Theoretical estimates, on the other hand, suggest that all non- precessional solutions should be safe from unpinning in the 
inner crust. Alternatively, one can note that an observation of a non-precessing star radiating gravitationally at fl and 2f2 
would imply that the bounds on the critical unpinning velocity from glitches and precession (equations [90] and I92[) are not 
reliable. 



5.2 Changes in shape 

Now let us examine the change in shape of the star in response to the pinning forces. We will consider the cases of elastic and 
magnetic deformation separately. 

Let us begin with a star deformed by a magnetic field of strength ~ B. The deformation A7 will be sourced by a 
magnetic stress of order B 2 , and, to a rough appr oximation, t he fractional deformation A// 1 will be of order of the ratio of 
magnetostatic and gravitational binding energies (| Jonesl l2002f ): 
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M B 2 R S 

I ~ GM^/R' ( ' 

The Magnus force will perturb the star, resulting in a magnetic stress perturbation of order BSB, and so a magnetic force 
per unit volume ~ BSB /R. This balances the Magnus force per unit volume of p SF RQ. 2 AI/I SF . Combining these results, and 
simplifying using 7 SF ~ p SF R 5 leads to 

sb m 2 

B ~ GM^/R' ( ' 

This is the ratio of the rotational kinetic energy to the gravitational binding energy, and will be much less than unity for all 
but the most rapidly rotating stars. We can therefore conclude that, in the case of a magnetically deformed star, the change 
in magnetic field, and therefore shape, caused by the pinning forces is small. 

A very similar calculation applies is we assume the pinning occurs in a crustal shell of thickness AR, shear modulus p, 
subject to a strain u. The deformation AI/I will be sourced by an elastic stress of order pu, and, to a rough approxim ation, 
AI/I will be of order of the ratio of Coulomb binding energy to gravitational binding energy, multiplied by the strain ( JonesI 
20021) : 



^L^ ^ AR U (97 ) 
I GM 2 /R (y ° 

The Magnus force will perturb the crust, resulting in an elastic stress of order pSu, and so an elastic force per unit volume 
of pSu/ AR. This balances the Magnus force per unit volume p SF RD. 2 AI/I SF . Combining these results, and simplifying using 
I SF ~ p aF R 4 AR leads to 

Su AR in 2 

u ~ R GM 2 /R' { ' 

This is a factor AR/R smaller than the corresponding magnetic result, showing that, in the elastic case, the pinning forces 
have minimal impact on the strain distribution. 

We can therefore conclude that pinning forces do not significantly alter the shape of our star. This means there is no 
danger of the pinning inducing fracture, and also that the gravitational wave field calculated previously is not significantly 
modified by Magnus force-induced shape changes. 



5.3 Superfluid instabilities 



We can also ask whether or not a superfluid vortex instability might be induced, of the form investigated recently by several 
authors ( Peralta et al. 20061 : Siderv. Andersson fc Comer 2008). The instability occurs in stars consisting of two fluid compo- 
nents: a superfluid neutron component, and a charged component, coupled via the mutual friction force. When the coupling 
between superfluid vortices and the star's charged components is su fficiently strong, a relative flow of super f luid a long the 
vortex array induces an inertial plane wave instability. In particular, iGlampedakis. Andersson fc JonesI l|2008i l2009h showed 
that such an instability would be triggered in a precessing star, if the wobble angle were sufficiently large. Could such an 
instability be triggered by the angular velocity mismatch in our non-precessing star? 

Taking the strong pinning limit of Glampedakis. Andersson fc Jones ( 20091 ) (their equations (46) and (47)) and writing 
the relative flow of superfluid along vortices, w, as w w £IRAI/Isf as given in equation (|85|l . we find that there do indeed 
exist unstable plane wave solutions; see appendix [B] for details. These occur for wavelengths smaller than a critical value: 
A crit w itRAI/Isf: 



\ critical _ on 

A ~ 30 cm 



R 



10 6 cnJ ^ 10- 7 J \I SF /I J 



This wavelength is longer than the inter-vortex separation in all stars of interest ijSaulsll 19891 ): 

1/2 



Ax v 



10" 2 cm — 



(99) 



(100) 



leaving a range of wavelengths where the instability could be operative. As shown in appendix [Bj the fastest growing pertur- 
bation has a wavelength of approximately A crlt /2, and a corresponding growth time of 

o is(— — ] ( ] !^)j^=n! | fi^rr) • inn) 



/ 



T 



0.1 



R 



io- 



In this equation x v denote the ratio of the proton to neutron densities, and the temperature T enters as the instability is 
sensitive to the temperature dependent electron-electron shear viscosity. These results indicate that the superfluid instability 
may be operative in stars of gravitational wave interest. If so, the instability could potentially unpin the vortices and, 
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presumably, eliminate the misalignment between superfluid and crustal principal axis necessary for our wave generation 
mechanism, e liminating the multiple har monic gravitational wave emission. 

However, van Hoven fc Levin (2008) have recently argued that the inclusion of hydromagnetic forces acts to stabilise the 
instability. They found that, when such forces were included, the relative velocity needed for instability was (BB CT /np SF ) 1 / 2 , 
where _B cr ~ 10 15 G is the critical magnetic field confined to the flux tubes. Combing with our relative velocity w m QRAI / I SF 
allows us to estimate the spin frequency above which this inertial wave instability would operate: 

/ =i 3 x 10 4 Hz («~) (^) (Qf) (102) 



R J \AI/I I \ 10 



This is very high, indicating that the magnetohydrodynamic forces might indeed be effective in suppressing the instability in 
all but the most weakly magnetised stars. 

To sum up, the multipole moments Q21 and Q22 that g enerate our gravitational wave fi e ld wil l be bounded by the crust's 
finite shear modulus and breaking strain, as described by Ushomirskv. Cutler fc Bildsten ( 2000T) . There are other possible 
failure mechanisms connected with the finite strength of the pinning and a possible vortex instability. Depending upon the 
detailed (and poorly understood) physics of the amount and location of the pinned superfluid and the pinning strength, these 
failure mechanism could potentially prevent the creation of a steadily rotating Q21 mass quadrupole. Conversely, observation 
of gravitational radiation at the two frequencies (fi, 20.) would imply that these failure mechanism are not operative. 



6 OTHER FREQUENCY SPLITTING MECHANISMS 

There exists (at least) two other ways in which an apparently steadily spinning star might emit at both Q. and 2fi, even 
without a pinned superfluid component. We will comment on these mechanisms briefly, and argue that they are much less 
efficient than superfluid pinning in producing the double frequency emission. 



6.1 Free precession of a biaxial star 

Consider a perfectly biaxial star, i.e. one with moment of inertia tensor of the form 

lab = IaSab + AlUaUb- 



(103) 



As is well known the free precession of such an object is very simple, consisting of the superposition of two steady rotations: the 
symmetry axis n a rotates about the fixed angular momentum vector with period ~ P, with a superposed slow (Pf p ~ PIo/AI) 
rotation about the symmetry axis ( Landau fc Lifshitz 19761) . If the pulsar beam lies exactly along n a , and if the beam is 
exactly axisymmetric about n„, then the radio pulsations will be steady at period P, with no processional modulation (Jones 



& Andersson feoOll) . The gravitational w ave emission is insensitive to the superimposed slow rotation, an d has the form given 
by equations I|50[) - ll53p. with J31 = A7 JZimm erman n fc Szedenits. Jr". Ill979t I Jones fc Anderss"onll2002l). This is the sort of 



motion modelled in 



Gal'tsov. Tsvetkov fc Tsirulevl (|1984l ). iBonazzola fc Gourgoulhonl i|l996h , iMelatos fc Pavnd (|2005f ). and 



Vigelius fc Melatosl (2009), where a magnetic field sources the deformation AI; see section [7] for discussion 



However, this is clearly a very special case. Firstly, it is unlikely that the moment of inertia tensor is biaxial around 
the magnetic axis. Crustal contributions are likely to be misaligned fr om the magnetic ones, particular ly if the crust retains 



any memory of a relaxed shape from a previous higher rotation rate (|Culter. Ushomirskv fc Linkll2003f ). and any significant 
contribution to I a b not of the form AIn a nb would make the inertia tensor triaxi al, breaking the symmetry and introducing 
detec table periodicities into the observed rotation (jZimmermann fc Szedenits. Jr ] 1 19791 : IZimmerman n 1980; Van Den Broeckl 
20051) . Secondly, pulsation profiles arc known to be non-axisvmmetric, as evidenced bv the complexity of pulse profiles (Lyne 
fc Graham-Smith 1 19981). Finally, the precession motio n would damp on some dissipative timescale, so would not be long-lived 
( Sedrakian. Wasserman fc CordesHl999t ICutlcr 2002). We therefore feel that such a motion is unlikely to be common in the 
pulsar population. 



6.2 Electromagnetic torques 

Magnetised neutron stars are acted upon by an electromagnetic torque. The torque consists of the familiar spin-do wn torque, 
scalin g as fi 3 , and also a so-called anomalous torque, a factor c/ (Q.R) larger than the spin-down torque, scaling as Q 2 (jGoldreich 
1970h : 



T a — j^j(ftd m d)e a bcOb m c 



(104) 



where m a is the star's dipole moment and a a factor of order unity, which Goldreich ( 19701 ) sets to unity, but Good fc N j 
1 19851) correct to —1/5. This torque is orthogonal to f2„ and so does no work on the star. However, it can have an affect upon 
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the rotation (|Goldreichlll97ol ; llVIelatosI I2OO0I) . To see this, look at the special case of a star with a biaxial deformation AI 
caused by crustal strain and a biaxial deformation AJ B caused by magnetic fields: 



hb = loSab + A/ c nC n c + A/ c n B n B , 



(105) 

where the unit vectors and n B define the symmetry axes of the crustal and magnetic deformations, respectively. The Euler 
equation of motion in the presence of a torque T a is 



dJa O o 7 

at 



T a . 



Combining the above equations and selecting the non-precessional solution by setting dj a /dt = we obtain 



AI c (n^^)e abc ^n^ + AI C (n^^)e abc ^n^ 
where m a = mnj and 

2 



A{il d n d )e abc il b n c , 



A = 



am 



As noted bv lGoldreichl l|l970l) . this torque can be absorbed into the moment of inertia tensor: 



e ab ^[AI c {n^)n^ + (A/ B - A)(n B «£)n B ] = 0. 



(106) 



(107) 



(108) 



(109) 



(110) 



It follows that the non-precessional solution is misaligned with the principal axes of I a b by an angle of order 
50- 

This misalignment will create a non-axisymmet ric mass distribution and lead to g r avita tional wave emission. That electro- 
magnetic torques have this effect was notes by ISedrakian. Wasserman fc Cordesl l|1999l) , who considered the effect of the 
non-anomalous spin-down torque; the principle is exactly the same. 

To gauge the significance of this, note that the magnetic moment m is related to the (surface) magnetic field strength 



via m ~ BR 3 , and, the magnetic deformation scales as A/ B ~ B 2 R 6 /(GM), so that 



.4 , 



R 



(111) 



(112) 



where -Rsch ~ GM/ (Rc ) is of order the star's Schwartschild radius. If follows that the misalignment angle is of order 

, A/ B R sdb 
AI R ■ 

It follows that this misalignment is very small unless the magnetic field makes a significant contribution to the moment of 
inertia tensor, and even in the case A7 B ~ A7 C ~ A7 the angle is suppressed by the factor Rscti/R ~ 0.2. We can therefore 
conclude that the anomalous torque is not as effective as pinned superfluidity in modifying the gravitational wave spectrum. 

In fact, there is a further caveat to be attached to this. The above argument shows that there is a misalignment between 
the spin and principal axes in a steadily spinning star. If the gravitational wave emission is then calculated by integration over 
the mass distribution, emission at frequencies Q and 2Q would be found in the standard way. However, one should presumably 
include the contributions to the stress-energy tensor from the magnetic field itself in this calculation. One might worry that 
inclusion of this additional component might eliminate the fi harmonic. This is clearly an interesting issue, which we will 
address in a future study. 



7 DISCUSSION 

Most targeted gravitational wave searches to date have been motivated by the possibility of the stars having 'mountains', 
i.e. non-zero values of the Q22 mass multipole moment. Such a deformed star, rotating rigidly about a principal axis, would 
emit at exactly twice the spin frequency. Radiation at frequency Q has received relatively little attention. In large part, this 
is because, in a single component rigid star, the necessary Q21 mass multipole moment would require the star to precess, 
modulating the observed timing of the radio pulses, and there has been little clear sign of such modulation in the radio 
pulsations of gravitational wave candidates. 

In this paper we have argued that a neutron star containing a pinned superfluid component can rotate in a steady 
non-precessional way, radiating gravitationally at both Q and 2f2, without there being any modulation in the observed radio 
pulsations. Fundamentally, this is because the pinned superfluid effectively acts as a gyroscope, 'sewn' into the crust of the 
star. This adds an extra piece to the angular momentum, allowing the system's total angular moment vector and rotation axis 
to coincide, even though no one of the crustal principal axes lies along the rotation axis. The superfluid pinning is crucial, 
and the axis of pinning must not be aligned with a principal axis of the moment of inertia tensor; if it is, the non-precessional 



© 2009 RAS, MNRAS 000.rflfT9l 



GW emission from superfluid neutron stars 15 



motion is about this axis, and only 20 radiation is produced. As we have noted, there exist a number of mechanisms which 
could break the pinning, depending upon the exact nature of how the superfluid vortices interact with the rest of the star. 

What does this mean for gravitational wave search strategies? Clearly, this motivates carrying out gravitational wave 
searches where data from two harmon i cs are combined to give a single detection statistic, as described for biaxial precessing 
stars by Jaranowski. Krolak fc Schut j ( 19981 ). The advantage of such a search over a single frequency search is that it would 



capture the full radiated signal. However, there is also a disadvantage: the two-harmonic signal is specified by 10 parameters, 
whereas the conventional single frequency search has only 7 parameters. It is therefore likely that the signal-to-noise threshold 
necessary to claim a detection would be larger in the two-harmonic search, as the larger parameter space increases the chance 
of a false detection. Ultimately, just how useful the two-frequency search is compared to the single frequency one will depend 
upon what sources nature chooses to provide. For a star that radiates only at twice the rotation rate, the two frequency search 
simply degrades our sensitivity to detection, without adding anything to the signal-to-noise. On the other hand, for a star that 
radiates appreciably at both harmonics, the two frequency search would surely be an advantage. According to the analysis 
of section T3.3I the two multipole moments contribute equally to the signal-to-noise ratio when IQ21I ~ 4| Q22 1 (neglecting 
the variation in the detector noise with frequency). In the extreme case of a targeted search, if the star radiates only via 
the Q21 multipole, the two frequency search would be essential. Clearly, some numerical experimentation to quantify when 
you gain and when you loose would be useful. One simple possibility, motivated by the less efficient nature of gravitational 
wave emission from the Q21 moment, would be to search at a single frequency, and only in the event of a detection look for 
harmonic structure. 

However, all of this begs the question: Why should the Q21 mass distribution be non-zero? What astrophysical scenario 
might produce such a state? Unfortunately, this is a far more difficult problem to address. For stars spinning rapidly enough to 
be of gravitational wave interest, the dominant force deforming neutron stars away from spherical symmetry is rotation, which 
is necessarily axisymmetric and can only contribute to the (non-radiative) Q20 mass multipole. It is therefore not obvious 
what might lead to a significant non-zero value of Q2i- However, exactly the same problem applies in arguing for the existence 
of a Q22 asymmetry. So, by this measure at least, there is no good argument to strongly favour either one of the two sorts of 

gravitational wave producing asymmetries over and above the other. 

One possib ility for producing asymmetry is a strong internal magnetic field dCal'tsov. Tsvetkov fc Tsirulevl 19841: Bonazzola 



Gourgourhon ll99ci lHaskell et allboosl; Iciolfi et alj|2009l: lLander fc Jones 2009). For instance, a magnetic field, symmetric 



about an axis along inclined at an angle x to f fle rotation axis, would produce a Q20 mass multipole when referred to a 
coordinate system with z-axis along direction , which, for < \ < 7r /2, translates into an ^-dependent linear combination 
of Q20, Q21, Q22 multipole moments, as given by the transformation properties of spherical harmonics under rotation: 

3 3 1 

Y20 = -p sin 2 x(^22 + Y2-2) pisinxcosxC^2i + Y2-1) + -(2 cos 2 x - sin 2 x) Y 20- (113) 

4v6 V6 2 

However, as is well known, the strong field neutron stars (the magnetars) all rotate too slowly to be of gravitational wave 
interest, while the more rapidly spinning stars (the young pulsars, the low-mass X-ray binaries (LMXBs), and the millisecond 
pulsars) have relatively weak external magnetic fields. Clearly, only if the internal magnetic field far exceeds the external one 
would magnetic deformations be of interest. 

Most of the known pulsars do indeed seem to be 'inclined rotators', i.e. have magnetic fields that are neither aligned with 
(X = 0) nor orthogonal to (x ~ i"/2) the rotation axis. What does this mean for their gravitational wave emission? The answer 
to this depends upon whether or not there is a pinned superfluid component. In the absence of such a component, there are 
two possibilities. Firstly, if the magnetic deformation were the only deformation, the star would precess, and either align or 



go or thogonal on some short timescale, leading to no gravitational wave emission or emission at only 2fl, respectively (|Cutler 



200 



3). Secondly, and more plausibly, there would be a corresponding crustal contribution to the inertia tensor, such that the 



in clined rotato r simp ly rotates about a principal axis of the total inertia tensor, radiating at 2f2 only. (As discussed in detail 
by Wasserman 1 20031 ). for a given inclination angle x an< A a given magnetic deformation, this would require a sufficiently large 



crustal deformation). However, if there is a pinned superfluid component, then the analysis of this paper applies, the motion 
is non-precessional, and emission at both il and 2f2 is possible, depending upon the orientation of the pinned superfluid with 
respect to the principal axes of the total inertia tensor. 

One particular scenario where magnetic fields might prove to be important is in accreting systems, if magnetic burial of 
the field takes place, resulting in stars with f ast s pins, strong internal fields and weak external ones. Such a scenario has been 
modelled in detail bv lMelatos fc Pavnel |2005l) and lVigelius fc Melatosl |2009l) , who showed that accretion of AM ~ 2 x 10" 3 M o 



was sufficient to generate a mass quadrupole of size AI/I » 2 x 10 -5 , biaxial about the magnetic axis. As noted by the above 
authors, this would generate a precessing star, radiating gravitational waves at the fi, 2f2 harmonics. 

However, in the absence of superfluid pinning, there is an objection to this scenario. The precessional motion would damp 
on a timescale much shorter than the timescale on which the magnetic mountain is built: the quadrupole is generated on a 
timescale related to the accretion rate M, roughly AM/M ~ 10 6 years for a typical LMXB, while the precession period, in the 
absence of pinned superfluidity, is of order P/e ~ 10 2 seconds. It follows that the quality factor of the precession would have 
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to be unrealistical ly large to allow the precession, and therefore the radiation at frequency Q, to survive (jJones fc Andersson 



200ll: ICutlerll2002h . This changes completely if a pinned superfluid is added to the model. In this case, the rotation remains 



close to the superfluid pinning axis, and there is no precession, allowing the inclined biaxial magnetic mountain to radiate 
continuously at the two frequencies. So, unless there exists some mechanism to maintain the precession, the superfluid pinning 
model of this paper is essential to realise the long-lived multiple frequency gravitational wave generation described bv Melatos 
& Pavne l|2005h and lVigelius fc Melatosl (|2009h . 



The solid crust is the other possible source of the deformation. As was computed in detail bv lCulter. Ushomirskv fc Link 



(2003), the crust in a pulsar is likely to have a Q20 density perturbation sourced by strains generated by the crust retaining 



a memory of a relaxed state at a higher rotation rate. If such a star is kicked, precession occurs, generating the familiar 
gravitational wave field. However, what would generate non-zero Q21 (or, for that matter, Q22) for our non-precessing stars? 
One possibility is that some of the glitches in young stars correspond to starquakes, i.e. sudden cracking events, and that 

these don't occur in completely axisymmetric way. 

In fact, a model of symmetry breaking in starquakes has been developed bv lFranco. Link fc Epstein (2000), who were 



attempting to explain the permanent increases in spin down rate seen after glitches in the Crab pulsar. Their suggestion was 
that the crust cracks in a non-axisymmetric way, so that, immediately after the glitch, the spin and principal axes are no 
longer coincident. They argue that a brief period of precession would follow, which, when damped, would leave the star with 
a permanently increased angle between the magnetic and spin axes. Crucially, to account for the systematic increase of this 
angle, Franco et al. argued that the magnetic field itself breaks the axisymmetry in the strain build up, and favours quakes 
along fault lines which tend to increase the spin axis-magnetic dipole axis angle. Clearly, such non-axisymmetric faulting 
could, in our pinned superfluid scenario, contribute to building non-zero values of the Q21 and Q22 multipoles. However, for 
this to work, the starquake would have to occur without significantly unpinning the superfluid, or at least the pinning would 
have to reform on a timescale short compared to the timescale on which the precession of the star is damped, or else the 
system would relax to rotation about a principal axis. 

The neutron star crust tecton ics model of Ruderman and collaborators provides another possible mechanism of symmetry 
breaking (|Rudermanll 197rl 1 1991allbh . This mechanism makes use of both the magnetic field and the core superfluid. The idea is 



that as a star spins down and its core superfluid vortices migrate outwards, they interact strongly with the magnetic flux tubes 
which would thread a type II superconducting core. This places a large strain on the crust. The exact deformation produced 
(i.e. the sizes of Q21 and Q22) would r equire a detailed calcul ation, and would depend upon the poorly understood interaction 
between the vortices and flux tubes ( Siderv fc Alpai 2009h . but Ruderman estimates that the corresponding stresses are 



large-large enough to shear the crust, and so could be a source of non-axisymmetry. 

To sum up, we have argued that pinned superfluidity allows a neutron star of spin frequency Q to emit steady gravitational 
wave signals at both Q and 2H, without there being any stellar precession. This motivates gravitational wave searches at both 
harmonics, even when targeting known pulsars with very smooth timing profiles. The non-detection of such radiation would 
imply that the necessary gravitational wave generating quadrupole moments aren't created in real stars, and/or the pinning 
isn't strong enough to support them. More interestingly, the observation of radiation at both harmonics from a non-precessing 
star would provide evidence in favour of pinned superfluidity within the star. 
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APPENDIX A: THE COUPLING TORQUE 



In our model, the superfluid component does not rotate steadily — its angular velocity vector Q^ F rotates about the fixed 
crustal angular velocity vector Q^. In this appendix we will show that the torque required to support this motion is supplied 
by the Magnus force which couples the superfluid and crustal com ponents via the superfluid vortices. 
The superfluid is acted on by a Magnus force (|Shahamll 19771) 



fa 



SF / SF 

: p e abc uj b (v c 



(Al) 



where u„ = e a b c ^b x c is the superfluid velocity, = e a bc^b x c the crust velocity, and tu a = 2f2„ the superfluid vorticity. 
We want to see what sort of torque this produces on the crust: 



tabcXbfc dV. 



(A2) 
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Combining the above equations and simplifying we find 

T a = 2e abc I p SF x b x d dVfl s /(n s c F - (A3) 
Jv 

In spherical symmetry (which we assume for the superfluid) the integral simplifies: 

/ p SF x b x d dV = S bd / p SF x 2 dV, (A4) 
Jv Jv 

while the moment of inertia tensor is, by definition, 



lab = P (r S a b - x a x b ) dV, (A5) 

■J V 

where r 2 = x 2 + y 2 + z 2 . Picking out a single component: 

/L F = / P SF (V 2 + z 2 ) dV = 2 [ p SF x 2 dV, (A6) 
Jv J V 

and so 

p SF x b x d dV = 5 bd I SF /2, (A7) 

where I SF = if F = Iy F = if F is the moment of inertia of the spherical superfluid. Inserting this into (|A3[) and writing in 
terms of the superfluid angular momentum jf F = I SF 0. a gives 

T a — e abc Q, b jf F . (A8) 

This is exactly the torque required to spin a superfluid of angular momentum jf F at a rate , confirming that the Magnus 
force does indeed supply the necessary torque to provide a self-consistent solution. 



APPENDIX B: THE SUPERFLUID INSTABILITY 

The short wavelength superfluid instability is described in detail in iGlampedakis. Andersson fc Jones! l|2009h . The complex 
plane wave frequencies that apply in the limit of strong vortex drag are given by equations (46) and (47) of their paper. The 
unstable mode has frequency 

! 1/2 

(Bl) 



1 2 1 
+ kw H uk H 

2 Xrj 



— (vk 2 x v ) 2 + ixp(x p — 1)Q SF vk 2 + S1sf(1 + x p) 2 — 2fl SF kwx p 



where x p is the ratio of proton to neutron densities, k the wave number, and v the electron-electron shear viscosity. Taking 
the limit of small wavelength and retaining the dominant terms in the imaginary part of the mode frequency we find: 



2Q 



SF 



2Q b 



fc 2 



Inserting the relation w ~ _Rfl SF A7"// SF gives 
2fi 2 f 2 R AI~ 



9f(a) 



k 2 k I SF 



(B2) 



(B3) 



The mode is unstable for Q(cr) < 0, so the plane waves are unstable for wavelength less than a critical value at which Q(cr) = 0: 

critical = nR A/ (B4) 

Further analysis of (|B3|) shows that the fastest growing mode has a wavelength ^ crltlcal /2. Substituting this wavelength into 
(|B3[) and inverting gives the mode growth timescale: 



_ <ivx v 

T = n 2 R 2 y ai 

This gives a good estimate of the growth time of the most unstable mode. 



(B5) 
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